Dipole oscillations of a Fermi gas in a disordered trap: damping 
and localization 
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Abstract. - We theoretically study the dipole oscillations of an ideal Fermi gas in a disordered 
trap. We show that even weak disorder induces strong damping of the oscillations and we identify a 
metal-insulator crossover. For very weak disorder, we show that damping results from a dephasing 
effect related to weak random perturbations of the energy spectrum. For increasing disorder, we 
show that the Fermi gas crosses over to an insulating regime characterized by strong-damping due 
to the proliferation of localized states. 



Introduction. — Ultracold atoms in disordered po- 
tentials are currently attracting considerable interest. 
They offer unprecedented possibilities to revisit many 
open questions on disordered quantum systems with accu- 
rate experimental control of relevant parameters and orig- 
inal measurement techniques. The spectacular progress 
achieved in disordered Bose-Einstein condensates (BECs) 
have recently lead to the first direct observation of Ander- 
son localization of matter- waves [1,2], showing remarkable 
agreement with theoretical calculations [3,4]. A future 
challenge to ultracold gases is the production of disordered 
Fermi systems. Even better than Bose gases, they would 
mimic systems of direct relevance for condensed-matter 
physics, such as dirty superconductors [5] and granular 
metals [6]. Moreover, ultracold gases provide original in- 
sights on transport phenomena without direct counter- 
parts in traditional condensed-matter physics. For in- 
stance, ultracold atoms trapped in harmonic potentials 
may undergo dipole oscillations - i.e. oscillations of the 
center of mass (CM) - that can be observed for tens of 
periods. Since dipole oscillations are undamped in pure 
harmonic traps, irrespective to the nature of particles, 
temperature and interactions [7], damping results from 
the influence of external potentials only. Hence, dipole os- 
cillations have been used to characterize several properties 
of quantum gases [8,9]. These include bosonic quantum 
degeneracy [10], inhibition of transport of Bose gases in 
one-dimensional (ID) optical lattices [11], fcrmionic band 



insulators in periodic lattices [12, 13], localization of spin- 
polarized fermions in disorder- free aperiodic lattices [14], 
and interaction-controlled transport of fermions in optical 
lattices [15], just to mention a few. 
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Fig. 1: Dipole oscillations of a 3D elongated (A = 8) disordered 
Fermi gas. a) The oscillations are excited by shifting the trap 
from z = d to z = at time t — 0. The disorder is also switched 
on at t = 0. b-d) Dipole oscillations, averaged over 25 realiza- 
tions of the disorder, for a speckle potential with correlation 
length <7r = 0.08 l z and various amplitudes Vr. The solid black 
lines are exact numerical calculations of Zcwiit) (Eqs. id),©), 
the dashed red lines are obtained from perturbation theory 
(Eqs. H},0,(jSJ)). Horizontal dotted lines are Eq. (JSJ). Here, 
we have used d= 5l z , T = and fi = 200 huj z [20] . 
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In this work, we study the dipole oscillations of a spin- 
polarized Fermi gas in a harmonic trap combined with 
a ID disordered potential. Dipole oscillations arc in- 
duced by a sudden displacement of the trap center (see 
Fig. Ola)), a technique which is routinely used in exper- 
iments with ultracold atoms. We find strong damping 
even for weak disorder, and identify a metal-insulator 
crossover when the disorder is increased. We show that 
for very weak disorder, damping results from dephasing 
of different modes oscillating with frequencies randomly 
shifted around the harmonic oscillator frequency. This 
leads to weakly-damped and almost centered oscillations 
(see Figs. QJb) andQJc)), which are accurately described 
by perturbation theory. For increasing disorder, the Fermi 
gas crosses over to a strongly insulating regime, character- 
ized by strongly-damped oscillations with large offset (see 
Fig. HJd)). This insulating regime is due to the prolif- 
eration of single-particle localized states. So far, dipole 
oscillations in the presence of disorder have been studied 
only in BECs [16,17]. In this case, damping is due to 
dissipation induced by repulsive interactions [17], a pro- 
cess completely different from the one we identify here for 
non-interacting Fermi gases. 

Framework. — We consider a Fermi gas at thermal 
equilibrium with temperature T in a 3D axially-symmetric 
harmonic trap of frequencies uj z in the longitudinal di- 
rection z and cu±_ in the radial directions x, y (the as- 
pect ratio is A = lu±/lu z ). The trap is initially cen- 
tered at (x,y,z) = (0,0, d). Because of the Pauli exclu- 
sion principle, the fcrmions populate the single-particle 
eigenstates of the displaced trap, ^n^,n y ,n-(x,y, z) = 
0n x .n H (x, y)tjn} (z) , according to the Fermi-Dirac distri- 
bution. Here, \4>n T ,n ) is the eigenfunction of the radial 
harmonic trap associated with the eigenenergy frui±(n x + 
n y + 1). The function ipif\z) = <p n (z — d) with tp n (z) = 
(n\2 n l z y^)-^ 2 cxp(-z 2 /2l 2 z )H n {z/l z ) is the eigenfunc- 
tion of the ID longitudinal harmonic potential associated 
with eigenenergy e n = ftw z {n + 1/2), where H n (z) is the 
Hermite polynomial of index n and l z = \J h/mu> z the 
oscillator length in the z-direction. At time t = 0, the 
trap center is abruptly shifted to z = and a ID ho- 
mogeneous disordered potential, V(z) = Vrw(z/ctr), with 
average (V) = 0, amplitude Vr and correlation length or,, 
is switched on (see Fig. QJa)). In the absence of disorder, 
this process induces undamped dipole oscillations of the 
CM along z, ^cm(^) = dcos(uj z t) [7]. In the presence of 
disorder, we write the CM motion, averaged over different 
realizations of the disorder, as 

ZcM(t) = Z osc (t) + Zoc, (1) 

where Z osc (t) ~ T(t) cos(uj z t) is the oscillating part with 
r(t) an envelope function giving the damping of the oscil- 
lations and Zoo is the oscillation offset. In the remainder 
of the manuscript, we evaluate the quantities T(t) and Zoo 
and identify their physical origin. 



Using the density-matrix formalism, the CM motion 
reads Z CM {t) = Tr[ze- i ' fct /' 1 p eS e +i ™/ h ], where 

is the single-particle Hamiltonian along the z axis, and 
p e ff is an effective density matrix. The expression of p e ff 
is obtained by tracing out the radial degrees of freedom, 
which is possible because the 3D Hamiltonian is spatially 
separable and unchanged in the x, y directions at t = 0. 
We find 

/w = E /St( y ,A) i^)<y£°i' ( 3 ) 

n 

where the effective Fermi distribution, f™ s (T,[j,,\) = 
Enllo el^^+^ll-MlAB^+i ' includes the occupation 
numbers of the radial-trap levels [18] and N = 
J2 n fefi(T> Pi A) is the total number of Fermions. In a 3D 
elongated Fermi gas (A = uj±/ui z > 1), the ID limit is 
achieved when the population of the transverse excitation 
modes can be neglected. This occurs for hu>± — fi k^T 
when (i > and fi k^T (degenerate limit) or for 
fy-u± ^> k-gT when fj, < and \fi\ ^> k-^T (classical limit). 

To evaluate explicitly ZcM^t), we use the eigenfunctions 
{\ip n ), n £ N} and the associated eigenenergies {E n } of 
Hamiltonian TL. The 2Tcm(£) can be decomposed as in 
Eq. HD with 

Z osc (t)= £ e-^- E ^ h p n J(iP p \z\iP n ), (4) 

n,p;E n ^E p 

and 

Zoo = Yl PeS {^pl^n}, (5) 

where the matrix elements p™g = (V ; n|/Seff|V'j>} include all 
the dependency on the parameters of the initial Fermi gas 
(T, /x, A) and the trap displacement d. The disordered po- 
tential affects the energy spectrum (i.e. both E n and \tp n )) 
of the harmonic oscillator, which, from Eq. ([4J, appears to 
have a twofold effect on Z osc (t). First, it induces random, 
incommensurate, energy shifts which dephase the different 
oscillating components of the sum in Eq. As we will 
see, this is the main contribution to the damping effect 
for sufficiently weak disorder. Second, it induces random 
modifications to the harmonic oscillator eigenfunctions, 
which affect the terms {tp p \z\i/j n ) and p^s- According to 
Eq. ([5|), these modifications are also responsible for the 
existence of an oscillation offset, Zoo. 

Equations (|4]) and ([5]) are the basis of both our analyti- 
cal and numerical calculations. In the numerics, we use a 
speckle potential similar to that used in many experiments 
on disordered quantum gases [1, 19], for which the auto- 
correlation function reads C(Az) = (V(z + Az)V(z)) = 
Vj|c(Az/ctr) with c(u) = (sinw/w) 2 . Figures [T](b-d) show 
dipole oscillations of a 3D elongated (A = 8) Fermi gas 
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at zero temperature for various amplitudes Vr [20] . These 
results show increasing damping for increasing disorder, as 
expected. More precisely, we find a crossover from weak- 
damping (metal-like) regime for very weak disorder (see 
Fig.[TJb,c)) to strong- damping (insulator-like) regime char- 
acterized by a significant offset for stronger disorder (see 
Fig. |T[d)). As we will show, weak-damping results from 
disorder-induced weak perturbations of the energy spec- 
trum, while strong-damping signals strong localization of 
the Fermi gas. 

Weak-damping regime. — We first evaluate analyt- 
ically Z osc from Eq. (fj|. We take into account first order 
disorder-induced shifts to the harmonic oscillator eigenen- 
ergies and approximate the eigenfunctions \ip n ) with the 
harmonic oscillator ones \<p^) . Note that, as we will show, 
the offset Z M requires second-order perturbation calcula- 
tions. For consistency in Eq. (JTJ, one should in princi- 
ple use also second-order perturbation for calculating the 
envelope Z osc (t). However, it leads to small corrections 
in Z osc (t) so we disregard them. We have checked nu- 
merically that perturbation of the eigenfunctions has a 
negligible contribution to the damping, at least until al- 
most complete dephasing has occurred. Then, Eq. ^ is 
significantly simplified since (ip p \z\ip n ) ~ (tp^ |i|</4°' ) ) = 
^{V n + l<5p,n+i + \fn5p >n -i). Using first-order pertur- 
bation theory on the eigenenergies E n only, we find 

Z osc {t) = l z V n/2(»+1)F^ +1 cos (w,t + 5V n t/h) , (6) 



where 8V n 



ordered potential, 5V n is a random quantity with (SV n ) 
and (6V*) = V^R n (a R /l z ), where 



Tt + l/ 



(c^VWI^) and 



For different realizations of the dis- 
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is related to the Fourier transform of the reduced correla- 
tion function of the disorder, c{q) = J -j=c{u)e~ lqu , and 



L%(x) are Laguerre polynomials [21]. When averaging 
Eq. ([6]) over realizations of the disorder, we use a Gaussian 
distribution in the variable 8V n with mean square fluctu- 
ation (SV 2 ) [22]. We then find Z osc (t) ~ dT(t) cos(u z t), 
where 



r(t) 



n>0 



'l)^S+l eX P 



Vjt 2 OR 

2h 2 L 



R, 



/ 0~R 

\l, 



(7) 

The envelope function T(t) is non-universal in the sense 
that its general shape does not depend only on the dis- 
order (i.e. on Vr, or an d t ne model of disorder, v) but 
also on the initial density matrix, i.e. on temperature (T), 
number of fermions (or equivalently, chemical potential /z) , 
trap geometry (A) and initial trap displacement (d). For 
instance, T(t) is in general neither an exponential nor a 
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Fig. 2: Damping time, r, versus Vr in ID and 3D (A = 8) 
Fermi gases, both for /i = 200 hu) z and <tr = 0.08 l z . The 
inset is r as a function of or for the 3D gas. The lines are the 
predictions of Eq. ((7} and the points are results of numerical 
simulations. The dashed line in the inset is the white noise 
limit. The parameters are d — 5 l z and T = 0. 



Gaussian function [23]. To compare Eq. (|7|) with numer- 
ical data, we introduce the damping time r, defined by 
the equation T(r) = 1/2. Although r depends on both 
disorder and initial state of the Fermi gas, one can obtain 
some universal properties from the fact that the contri- 
bution of disorder (which appears only in the coefficient 
of the exponential function in Eq. ([7])) is separated from 
that of the initial state (in the quantities Fjf^ +1 only). 
In particular, we find r cx ?i/|Vr| for any model of disor- 
der and initial density matrix. Numerical results for ID 
and 3D Fermi gases at zero temperature shown in Fig. [2] 
confirm this prediction. The behavior of r versus the cor- 
relation length ctr, is more complicated. In the white-noise 
limit, ctr -c Iz/^/np, where rf = fi/hoj z is the longitudinal 
Fermi level, there are many disorder peaks within the typ- 
ical wavelength l z / \fnp of the Fermi gas. The disordered 
potential felt by the wavefunctions lyn "*) with n < n-p al- 
most averages out. When or increases, the eigenfunctions 
are more sensitive to the disorder so that we expect that 
r decreases. This behavior is confirmed by Eq. (J7J): In the 
white noise limit the functions R n are independent of <tr 
so that t cx y/lz/crR- The coupling of the Fermi gas to the 
disorder is maximum when ctr ~ l z / ^Jtvf [24]. Beyond, we 
thus expect that r increases. In the limit of large correla- 
tion length, or ^> l z y/2np, the typical size of the disorder 
peaks exceed the size of the Fermi gas. Hence, the disor- 
dered potential induces only an energy shift which is ap- 
proximately equal for all the cigenstatcs below the Fermi 
energy, and we recover harmonic oscillations (i.e. r — > oo). 
This is again confirmed by Eq. ([7]): Since the width of c 
is unity we find R n {o- K /l z ) cx {l z /cr R ) 5 and r cx (ctr// z ) 2 . 
These features explain the non-monotonous behavior of r 
versus or obtained in the Inset of Fig. [2l 

In contrast to damping, the oscillation offset results 
only from perturbation on the eigenfunctions \tp n ) (see 
Eq. ([5])) [25]. Since the average of the disorder vanishes 
everywhere, (V(z)) = 0, first-order perturbation is not 
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Fig. 3: Oscillation offset, Zoo, versus Vr (same parameters as 
in Fig. [2}. The points are results of numerical simulations and 
the lines are Eq. (JSJ). The inset is Z x as a function of <tr for 
the 3D gas. The dashed line is the white noise approximation. 

sufficient to calculate the offset. We thus evaluate Eq. ([5]) 
up to second order, and find 

n>m 

with 

^(r)^/^< v5 r K ) < " i) """ +li "-' ( ' ;2 » >< 

e-«" [Li(K 2 ) - + - ii-,(« 2 )] ■ 

As for damping, the behavior of the offset as a function of 
Vr is simple: we find oc (Vr/H(jJ z ) 2 for any initial state 
of the Fermi gas and trap displacement. This is confirmed, 
for small enough values of Vr, by the results of numerical 
calculations plotted in Fig. [3] Again, the dependence of 
Zee versus <tr is non- monotonous, as shown in the Inset 
of Fig. [3J In the white-noise limit, <jr <C l z /y/nF, we 
find that Z^ oc ctr, while for large correlation length, 
cr 3> l z \ / 2np, we have Zoo oc I z /<Jr. 

Let us comment on two important properties. First, we 
find that strong damping occurs even for very weak disor- 
der. For instance, Vr = hfrio z (equal to 0.025/x in Figs. [3J 
and 2]) corresponds to r ~ 2tt/uj z , meaning that damping 
occurs typically on few oscillations (see Fig. U|c)). Sec- 
ond, we find that, for the same chemical potential, the 
effect of disorder is weaker in ID than in 3D (the damping 
time r is larger and the offset Z M smaller). This can be 
understood intuitively from the properties of the effective 
Fermi distribution f™ s (T, /x, A). Indeed, for stronger ra- 
dial trapping (i.e. for increasing A or in ID compared to 
3D) and fixed Fermi energy, the population of low-energy 
states, which are more affected by the disorder, decreases. 
Similarly, we have found that finite temperature increases 
the damping time r and decreases the offset Z^ although 
the effect is quite weak. For instance, we find that a tem- 
perature T = 0.5 Tp, where Tp is the Fermi temperature, 
leads to a small correction of about 10%, without affecting 
the general behavior of r and Z^ versus the parameters. 



Note also that the assumptions at the basis of our per- 
turbation approach are justified, a posteriori, by the ex- 
cellent agreement between numerical calculations (points 
with error bars) and perturbation theory (lines) in Figs. [5] 
and [3J This is even more striking in Fig. [TJb,c) where 
we plot the dipole oscillations as obtained indepen- 
dently (i) from exact numerical calculations (based on 
Eqs. (H]), ([5]); solid black lines) and (ii) from the predic- 
tion of Eqs. IP),©,© (dashed red lines). They are hardly 
distinguished in Fig. [TJb,c), indicating that perturbation 
theory is indeed very accurate as long as oscillations are 
visible (i.e. for Vr < 5huj z for the parameters of Figs.[T][3|). 

Strong-damping regime and localization. In 

the strong-damping regime (i.e. for Vr > 10 hui z for the 
parameters of Figs. [T][3]), dipole oscillations arc largely 
suppressed (see Fig. QJd)) and it is difficult to define an 
envelope function T(t). At the same time, the offset be- 
comes significant (see Fig. [3]) and tends to ~ d for 
very large Vr [26]. This behavior signals the crossover 
to a strongly insulating regime which can be related to 
the onset of single-particle localization. Indeed, as shown 
by Eq. ([5]), the oscillation offset Z M is the sum of the 
average position of each eigenfunction in the disordered 
trap, z n = (ip n \z\ip n ), weighted by the corresponding pop- 
ulation, p™^™ [25]. Extended states -which are centered 
around the trap minimum (z n ~ 0)- do not contribute to 
Zqo. Conversely, localized states which are spread at ran- 
dom positions in the trap, z n ^ 0, may contribute to Z^, 
depending on their relative population 

The two important features that determine the oscil- 
lation offset are thus the localization properties of the 
eigenfunctions in the presence of the harmonic trap, and 
their populations, which are governed by the initial, dis- 
placed Fermi gas. In order to interpret the dipole motion 
in the strong-damping regime, we now discuss these fea- 
tures. The eigenfunctions \ip n ) arc obtained by numerical 
diagonalization of Hamiltonian TL and characterized by 
two quantities: (i) the localization center, z n , and (ii) the 
participation length, P n = 1/ J dz\ip n (z)\ 4 , which gives 
the typical extension (width) of the quantum states [27] . 
Their populations are obtained by projecting the ini- 
tial state of the Fermi gas on the eigenfunctions \ip n )- Fig- 
ure0]shows all these features for a ID Fermi gas at T — 0, 
with /i = 200Hlu and initially centred at d = 51 z . For 
each eigenenergy E n , we plot a vertical line of length P n , 
centred at position z n , and weighted by (color scale). 
The different Figs.[3Ja)-(f) correspond to the same realiza- 
tion of a speckle potential but different values of Vr. As 
can be anticipated, in the weak-damping regime, we find 
that the eigenfunctions are not localized (see Fig. IUa)). 
The values of P„ are close to the participation length P° 
calculated for the harmonic trap without disorder (solid 
blue line). The corresponding z n slightly fluctuate around 
the trap center, so that Z M <C d. Increasing the ampli- 
tude Vr of the disordered potential (see Figs. [4|b) and 
SJc)), most of the populated states remain extended but 
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Fig. 4: Localization in a disordered trap. Vertical lines centred in z n and of length given by the participation length P n are 
plotted at energy E n . The color scale is given by the weight p™g": it goes from zero (white regions) to one (dark regions). For 
clarity, the background (BG) is in light blue color. The horizontal dashed lines indicate the spatial region covered by the initial 
Fermi gas, i.e. d — l z s/1rw < z < d + l z s/2np- The different figures refer to different values of Vr and the same realization of 
speckle disorder. The solid blue line in Fig. (a) shows the participation length Pi°' of harmonic oscillator eigenfunctions in the 
absence of disorder (Vr = 0). 



become more and more perturbed. In the crossover regime 
(Figs Eld) and[3Jc)), strongly localized states, randomly 
distributed in the trap, appear at low energies. The states 
with lowest energy, which are more sensitive to the disor- 
der and thus more strongly localized, are populated when 
they lie in the spatial region covered by the initial Fermi 
gas (delimited by the horizontal dashed lines in Fig. [J). 
These states, being spatially separated from each other, do 
not contribute significantly to the sum that determines the 
oscillation (see Eq. ([4])). In contrast, they do contribute to 
the offset (see Eq. ((SJ) since they are mostly located on the 
side of the initial gas. The populated states with higher en- 
ergy are not localized and extend beyond the spatial region 
covered by the initial Fermi gas. They can thus contribute 
to the oscillation, but not significantly to the offset since 
their average center is close the trap center. Finally, for 
very large amplitudes Vr of the disorder (Figs. S{f)), the 
Fermi gas enters the strong-damping regime. In this case, 
most of the populated states are strongly localized and lo- 
cated in the region covered by the initial Fermi gas, with 
energies about up to the chemical potential. Note that the 
upper part of Figs. rj|e) and0|f), above the dashed line, 
shows states which are strongly localized but not popu- 
lated since they lie outside the initial Fermi gas. Some 
extended states with energy above the chemical potential 
can also, appear but with very small population. The CM 
of the Fermi gas is thus frozen at Z a fi ~ d. This explains 
qualitatively the behavior observed in Fig. [3l 



As we have shown, localization is a crucial ingredient for 
the metal-insulator crossover investigated here. Although 
localization is mainly due to the disordered potential, it 
is clear that the harmonic trap significantly affects the lo- 
calized states. For instance, by virtue of finite-size effects, 
the lowest-energy states are localized near the trap center 
(see Figs. fJJ. Moreover, the exponential decay of local- 
ized states, which is the most striking signature of An- 
derson localization [1,2], is strongly suppressed since the 
harmonic potential dominates at large distances. More 
surprisingly, our numerical data show that localized and 
extended states can coexist in the same region of the spec- 
trum at intermediate energies (see for instance Figs. HJd) 
and0Jc)), an effect due to the trapping potential. In the 
future, it will thus be interesting to study single-particle 
localization in disordered traps. On one hand, it differs 
qualitatively form localization in homogeneous disorder, 
and, on the other hand, it is directly relevant to experi- 
ments on disordered quantum gases. 

Conclusions. — Dipolc oscillations are an important 
tool for studying the dynamical properties of ultracold 
gases [8-17]. We have shown that dipole oscillations of 
trapped, disordered Fermi gases reveal a metal-insulator 
crossover for increasing disorder. Weak disorder induces 
weak-damping associated to weak perturbations of the en- 
ergy spectrum. Stronger disorder leads to strong- damping 
characterized by a large oscillation offset, which signals the 
onset of localization. We have related the insulating prop- 
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erty of the Fermi gas to localization of energy eigenfunc- 
tions in the disordered trap. We have provided analytical 
predictions in the weak-damping regime and numerical re- 
sults in the strong-damping regime using experimentally 
realistic parameters for both ID and 3D gases. Our pre- 
dictions can thus be readily tested experimentally with 
ultracold Fermi gases, in speckle potentials or other mod- 
els of disorder. Our results can be easily extended to finite 
temperature Fermi gases, but we have actually found weak 
differences compared to the zero-temperature case. 
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